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ABSTRACT 


The development and behaviour of bed forms in alluvial channels 
was studied by analyzing the results of a series of labératory flume 
vests’ 

A tilting flume abt a sand bed was employed for a total of thirty- 
one runs. An analysis was subsequently performed to determine the rela- 
tionships between the measured variables at conditions near the transition 
from one bed configuration to another. 

The analysis resulted in the representations of bed-configuration 
boundaries upon two graphical relationships. The first relationship 
involved a densimetric form of the Froude number, the concentration of the 
bed-material discharge and the ratio of the average depth of Flow to the 
median diameter of the bed-material particles. The second relationship 
was identical to the first with the exception of the Froude number which 
was replaced by the slope. 

The results were then compared with existing methods of predicting 


bed configuration. 
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SYMBOLS 


Channel width (ft.). 


Concentration of bed material discharge by weight of total 
discharge (parts per hundred thousand, ppht). 


Velocity of propagation of a bed wave (ft./sec.). 


Fall diameter of bed material (mm. or ft.). 
Median diameter of the bed material (mm. or ft.). 
Densimetric form of Froude Number FR' = pp /yth 


Constant of gravitational acceleration PEE 
Average depth of flow (ft.). 

Average discharge per foot of channel width (cfs/ft.). 
Discharge of water-sediment mixture (cfs). 


Slope of the energy-grade line equal to the slope of the 
water surface cr the bed slope during equilibrium conditions. 


Average velocity according to the continuity principle (ft./sec.). 


Variable describing physical characteristics of the bed material 
(gradation, particle shape, etc.). 


Dimensionless parameter V. = 3/vg 050 
v 


Buoyant unit weight of the bed material (1b./ft.). 
alah aati eo 

Height of bed wave (ft.). 

Length of bed wave (ft.). 

Kinematic viscosity of water-sediment mixture teva 


Density of bed material particles (slugs/ft.7) 


Density of the fluid belnasihtec): 
Shear stress acting on the bed RSI A aD 


Fall velocity of bed-material particles (ft./sec.). 
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CHAPTER | 


INTRODUCTION 


1-1 General 


Alluvial channels, which may be defined as being capable of form- 
ing their own boundaries by erosion, transport and deposition of sediment, 
have been used, modified and even created by man for centuries. However, 
with all his experience, man has not yet been able to inteesern the 
mechanism of these processes completely. 

Because the fluid in alluvial channels has two constituents (water 
-and sediment), the bed is characterized by various configurations which 
form according to the flow conditions. The determination of the rather 
ill-defined boundaries between different types of bed form is basic to 
the quantitative study of sediment transport and resistance to flow in 
alluvial channels. 

Studies in the specific area of bed forms can be classified under 
two headings. The first is observation of the phenomenon (both in the 
field and in the laboratory) and the subsequent derivation of empirical 
formulas. The second is derivation of formulas from assumptions or theories. 
The latter method seems to have been less successful in describing the 
overall phenomenon. 

Until recently, few experimenters noted the form of bed roughness 
at known flow conditions, so work under the first heading is relatively 
backward. Therefore, laboratory investigations are necessary to define, 


as completely as possible, the presently sketchy bed form phase boundaries. 
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12 Aim and Scope of Study 


The object of this study has been to define more precisely the 
conditions at which transitions in bed form configuration occur, and to 
determine the extent of any transitional areas which may exist. 

The experimental procedure consisted of the adjustment of any 
two flow parameters (primarily slope and depth) in such a manner as to 
approach the phase boundary with due attention to changes in bed form. 

The boundaries were approached from both sides to ensure the determination 
of transitional zones. Measurements were made of as many basic quantities 
as possible (including the bed form dimensions, velocity profiles, etc.). 
Although these measurements may not all be directly related to this part- 
icular thesis, they may provide data for later work without duplication of 


experimental efforts. 
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CHAPTER II 


LITERATURE REVIEW 


Probably the first attempt to determine an empirical formula to 
describe the movement of a bed form was that of Sainjon (ref. 4, p. 12). 


His studies on the Loire River in France led him to suggest that: 
2 
c = 0.00013 (ve elas (2-1) 


where c is the velocity of propagation of a dune (metres/sec) and Mes is 
the velocity of water on the surface of the river. It is significant 
that this formula excludes the dune height and water depth. 

In 1883 Darwin (ref. 6, p. 175) published the results of experiments 
on sand ripple peta due to oscillatory motion of water over the bed. 
He felt that ripple formation was due to the generation of vortices in a 
layer between two uniform streams (water and sand). However, the first 
documented flume study was performed by Deacon (ref. 4, p. 12) in 1892. 
This author related the formation of ripples to water velocity. The 
experiment extended to the plane bed stage, and excellent descriptions 
were made of the mechanism of ripple formation. The most extensive 
early flume tests were those of Gilbert and Murphy (ref. 6, p. 178) in 
1914, and related slope, discharge and channel shape with sediment transport. 
Formation of all bed forms, including the antidune, was daseribed: 

Exner (ref. 6, p. 176) in 1932, employed a more analytical approach 
to the problem. He assumed that the capacity of water to transport sediment 
depends on velocity; therefore, acceleration results in scour and deceler- 


ation results in deposition. By assuming a uniform velocity distribution 
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and a sinusoidal bottom shape of wavelength A and amplitude "'a'', Exner 


concluded that: 


Cc = qk 


r ———~5- = velocity of propagation of ripples, and (2-2) 
(y-n) 
= + - 
» tr ss ead ( zs aaa) = bed elevation, (2-3) 
(y-n) 


where q is the discharge per unit width. The height of the bed is n and 
the height of the water surface is y (above some arbitrary datum). Although 
his theory does not indicate how the sinusoidal form is generated from a 
plane bed, it does explain the transformation from sinusoid to an asymmetric 
wave with a gentle upstream slope. Figure 2-1 is an illustration of Exner's 
mathematically derived dune profile. Exner further considered water 
surface slope and friction and derived complex relations for bed-form 
shape and channel width. 

In 1934 Kramer (ref. 4, p. 16) employed flume tests to determine 
the range of conditions for which a certain sand was suitable for hydraulic 
model tests. Rather than relating bed formation directly to velocity, as 
was the custom previously, Kramer used tractive force per unit area on 
the bed (or boundary shear) as a parameter. He concluded that the range's 
lower limit (formation of ripples) is dependent on the tractive force and 
the size and gradation of the bed material; the upper limit (formation of 
dunes) is dependent, in addition, on the hydraulic slope. In the ensuing 
discussion on Kramer's paper, Straub (ref. 6, p. 178) in 1935, concluded 
that initial movement and development of bedforms is dependent on the 
Froude number, gradation, and tractive force. Straub's experiments indicated 


that the height and sharpness of bed forms reaches a maximum at the critical 
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energy of flow and decreases as the energy increases. He also noted that 
the length of ripples is much less sensitive to changes in flow conditions. 
Shields, in 1936, investigated the papeg alae of bed forms. The 
graphical representation of his work (ref. 6, p. 101) has become known: as 
the Shields diagram (see Figure 2-2). Shieldsconcluded that bed formation 


was mainly controlled by Te where ‘c is the critical drag, y' is the 


2 


submerged particle density and d is the grain diameter, and by d where 6 
6 


is the thickness of the laminar film. 
Goncharov (ref. 6, p. 179) defines bed formation in terms of some 


non-transporting velocity Ver) and a relative tp dbed Ns The 


velocity at which bed forms begin to develop is V' where: 


1/12 
yo = 25(K) (2-4 
VON Yo | 
The velocity at which bed forms disappear is V''' where: 
1/12 
tii ia 
NON k 


The velocity at which maximum development occurs is V'' where: 
VUNG me Vo te), 25°" 3 (2-6) 
Bakhmeteff (ref. 4, p. 15) involved the use of the ''kinetic flow 
factor'' (square of the Froude number) which resulted in expressions for 
the velocity of propagation of bed forms such as: 


c= 0.0188 ve - 0.0292 gd (2-7) 
ci BY, 42 
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Velikanov's work (ref. 6, p. 181) in 1936, was mainly concerned 
with the structure of turbulence in open channel flow. He proposed a 
mathematical model explaining that turbulence (in the form of large scale 
periodic eddies) could cause scour and deposition on the bed. Furthermore, 
he was able to show that the generation of a sinusoidal bed form from a 
plane bed (a basic assumption in Exner's development) could be determined | 
from turbulence considerations. Similar work was performed by Tison 
(ref. 6, p. 184) who, in addition, concluded that periodic ripples do not 
occur in the laminar flow of a highly viscous fluid. 

Ismail (ref. 4, p. 20) studied flow in a closed conduit with bed 
forms, using Reynold's number as a criterion for ripple formation. He 
concluded that some critical Re exists, below which part of the sediment 
moves as bed load, and above which all] sand moves in suspension (dunes are 
swept away). 

Anderson, Kondrat'ev and Kennedy (ref. 6, p. 184) all considered the 
problem in terms of potential flow. Kennedy derived a complicated analy- 
tical solution which contains an unknown dimensionless parameter that is 
dependent on depth and velocity of flow and the fluid and sediment proper- 
ties. As such, his approach adds little to the overall understanding of 
the problem. 

Most recent work e in terms of the Froude number (in some form), 
as illustrated by Figure 2-3 from the work of Simons and Richardson (7). 
Subsequent work by these authors (10) relates the phases of bed forms to 
stream power (t V) and median fall diameter of the bed material. This 
relationship is shown graphically in Figure 2-4. This analysis employed 


both flume and field data and was intended to facilitate the prediction 
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FIG. 2-3. SIMONS AND RICHARDSON CRITERION FOR 
PREDICTING BED CONFIGURATION 
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of the form of bed roughness, given slope, depth, velocity and fall 
diameter. 

Dimensional analysis led Liu (5) to the use of a Reynolds number 
as a parameter in defining the various phases of bed forms. A comparison 
of Liu's criterion for ripple formation, Shield's initiation of motion 
curve and the Simons and Richardson curves was published by the last 
authors (8) in 1961. It is presented graphically in Figure 2-5. After 
discussion of Liu's paper had concluded, the author was of the opinion 
that ripple formation was caused primarily by an interaction of turbulence 
and instability of the sand-water interface. However, Vanoni and Brooks 
(12) felt that the consideration of stratified flow was inappropriate as 
the sand layer of ''flow'' has many characteristics which differ from those 
of a very dense fluid. Many subsequent analyses have been made with refer- 
ence to turbulence, but due to their complex nature and inherent assump- 
tions, they remain rather academic. 

Cooper (3), in an analysis of world flume data, indicate the phase 
boundaries on general solution surfaces for coordinates of Fr', C, h/D, 
and $, C, h/D. This more general, dimensional approach is more likely 
to result in a basic understanding of the problem. Their study reveals 
areas in which flume studies have not provided sufficient information 
and in which future efforts should be concentrated. 

The problem of identifying and classifying the physical features 
of bed forms is not nearly as controversial a subject as is the prediction 
of the form of roughness for a given flow condition. The classification 
system proposed by Simons and Richardson (7) has been generally accepted. 


Yalin (13) suggested that ripples and dunes share a geometric resemblance 
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ioc): SIMONS, RICHARDSON AND ALBERTSON 
CRITERION FOR PREDICTING BED CONFIGURATION 
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but that ripple size depends on grain size (not flow depth) whereas dune 
size depends largely on the depth of flow (and is independent of grain 
diameter). Simons and Richardson also indicate that the length of these 
bed forms is much less sensitive to changes in flow condition than is the 


height. 
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CHAPTER III 


EXPERIMENTAL PROGRAM 


Shed Objective and Scope of Study 


As stated briefly in the first chapter, the objects of this study 
were to define the conditions at which transitions in bed configuration 
occur, and to determine the nature of any transitional configurations 
which may exist. 

The experimental approach used in pursuing these objectives was 
to slowly change the slope and depth (to maintain uniform flow conditions) 
while keeping discharge constant. According to existing knowledge, the 
sequence of bed configurations which could be expected for increasing 
slopes (and decreasing depths) was: (1) ripples; (2) dunes; (3) upper 
flow regime. A sequence of bed configurations in the reverse of that 
stated above could be expected to result from decreasing the slope (hence 
increasing the depth). Upon commencement of transition from one bed to 
another, the flow conditions were allowed to stabilize and measurements 
were taken. Attempts were made to develop bed configurations according 
to both the sequences mentioned above, so as to determine the existence 
of any transitional bed forms me to determine the possible effects of 
antecedent bed forms on the development of the new configuration. 

However, the physical nature of the laboratory apparatus imposed 
some limitations on the experimental range of slopes and depths which 
could be achieved. This became critical in the case of high slopes as 
discharge was limited to approximately one cfs, and at lower discharges 


meandering tended to occur within the flume's four foot width. Thus, a 
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lack of testing resulted in the upper regime area (specifically, in defining 


the upper boundary of the transition from dunes to upper regime). 


Sia Apparatus 


While the aforementioned limitations did exist, the experimental 
apparatus provided a sufficient range of adjustment to allow attainment 
of the basic objectives of this study. 

The flume (which is illustrated in Plates 3-1 and 3-2, and in the 
schematic diagram of Figure 3-1) is located in the Graduate Hydraulics 
Laboratory of the University of Alberta. It has a width of four feet 
and a length of one hundred feet (although the use of baffles and a wave 
suppressor reduced the effective length to eighty-six feet). The entire 
water-sediment Are was recirculated through one of two andes lines. 
Gate valves regulated the flow so that at higher discharges a four inch 
diameter return line was used, while at lower discharge a two inch diameter 
line was employed. This procedure kept the velocities in the return lines 
high enough to prevent the development of bed forms (due to deposition of 
suspended material) which could have reacted to any changes in the flow | 
conditions. As a further precaution against external disturbances to 
the system, a baffle was installed at the outlet of the return line to 
reduce the turbulence at the entrance to the flume and a wave suppressor 
was employed to reduce wave action at higher discharges. The material 
used to pack the baffle was aluminum shavings which effectively reduced 
turbulence and resisted corrosion. 

The water-sediment mixture was recirculated by means of a Wilfley 
sand pump which had a maximum discharge of approximately one cfs. It was 


capable of pumping a slurry of up to 80% sand, a feature which became of 
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PLATE 3-1. THE EXPERIMENTAL FLUME 


PLATE 3-2. THE EXPERIMENTAL FLUME 
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use when the flume was refilled with water (causing bed material to slump 
into the inlet to the pump). 

The discharge was determined by means of a Foxboro magnetic flow 
meter and chart recorder. The device had two calibrated channels, one 
of which read from.zero to one cfs, the other of which read from zero to 
twenty cfs. The flow meter (shown in Plate 3-3) was installed vertically 
to eliminate air pockets or sediment deposition near the electrodes. 

Plate 3-4 shows the manometer board which was employed to determine 
the water surface elevation above pressure taps located every ten feet 
along the centerline of the flume bed. Plastic tubing ran from each tap 
to the board, which had an infinitely adjustable slope. The system was 
flushed periodically Mi water from the city mains so as to purge the 
lines of any air bubbles which could affect the manometer readings. 

As the entire water-sediment mixture was recirculated, it was 
desirable to take a total load sample at the outlet of the return Jine. 
This was accomplished with the use of a slice sampler. This device was 
mounted on a threaded shaft which, when rotated, caused the sampler to 
traverse the flow as it exited the return line vertically downwards. It 
was designed to collect a sample of approximately five hundred ml (de- 
pending on the cranking speed and water velocity). An enlarged bottom 
portion of the sampler was designed to act as a reservoir to prevent 
overflow due to the high entry velocities. The width of the opening was 
made several times larger than the median grain diameter to prevent sedi- 
ment from becoming lodged in the entrance. 


The dimensions of bed forms and the slope of the bed relative to 
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PLATE 3-3. MAGNETIC FLOW METER AND CHART RECORDER 


PLATE 3-4. MANOMETER BOARD 
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19 
the flume bottom were determined by the bed plotter shown in Plates 3-5 
and 3-6. This device consisted of a probe which was controlled so that 
it remained at a constant distance from the bed while the carriage upon 
which it was mounted traversed the length of the flume on rails parallel 
to the flume bed. The head of the probe contained an active electrode, 
at the bottom tip, and a remote grounded electrode. When the tip of 
the probe approached a non-conductor (such as the sand bed) the impedance 
between the two electrodes increased. This change in impedance caused a 
servo motor to restore the original separation between the probe and the 
bed. The movements of the probe and the carriage were monitored by potentio- 
meters which caused the bed profile to be recorded on an X-Y plotter. 

A Sanborn velocity probe was also mounted on the carriage mentioned 
above. It was used to obtain velocity profiles at the centerline and 
third iars of a cross-section during the run. However, due to the nature 
of the probe, the lowest point cn the profile was 0.03 feet above the bed 
and the highest point was 0.03 feet below the water surface. For this 
reason velocity profiles were not recorded for runs with shallow depths 
(less than about 0.2 feet). 

The bed material, while actually being one of the variables of 
the study, remained the same for all the tests, so it will be considered 
as having been part of the apparatus. The bed was comprised of Selkirk 
sand with a specific gravity of 2.65 and median diameter of 0.26 mm. The 


grain size distribution of this material is shown in Figure 3-2. 


3=2 Procedure for Determining Slope and Depth 
Due to the existence of bed forms on the bed, it was necessary 


to develop a reliable means of measuring depth and slope. It was 
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PLATE 3-5. BED PLOTTER AND CARRIAGE 


PLATE 3-6. BED PLOTTER CONTROLS 
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impossible, because of the undulations in the bed, to measure h and § 

over a short reach while the run was in progress. Furthermore, in order 

to prevent the bed forms from being swept away at the termination of a run, 
the depth was slowly increased (the flume was gradually flooded), thus 
changing the quantity of water in the previously closed system. Therefore, 
extensive use was made of the bed plotter and manometer board to determine 
the quantities h and S, 

Figure 3-3 is a schematic diagram of the manometer board during a 
typical run. The water surface is denoted by line 1-1, which has a slope 
So. The line 2-2 represents the water Sit race after completion of the run 
(with the flume level). Line 2A-2A represents the surface after correction 
for the flume slope during the run. Finally, line 3-3 indicates the loca- 
tion of the bed, a distance h' below line 2A-2A. This distance h' was 
determined by measuring, with a point gauge, the depth of water over a 
known point, A, on the bed profile (which had previously beaudleve died off 
with a trowel). The X-Y plot of the bed profile (such as that shown in 
Figure 3-4) was used to correct the depth, h', relative to the mean bed 


elevation. The plot also served to correct the bed slope, S,, should it 


we 
differ from the slope of the flume. Thus, the depth during the run was 
known to be the difference between the manometer readings of line I-] 


(the water surface during the run) and line 3-3 (the bed location during 


the run). 


a= 3 Classification of Bed Configuration 


As the main objective of this study was to classify the various 
types of bed configuration, it was necessary to establish a reliable 


means of identifying each bed form. 
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FIG. 3-3. ILLUSTRATION OF USE OF MANOMETER 
BOARD FOR DETERMINATION OF DEPTH 
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This was primarily accomplished with a visual inspection of the 
bed, as described more fully, with the aid of photographs, in section 
4-2, Unfortunately, the configurations are not as readily identifiable 
in the photographs as they were in the flume. 

The visual identification was largely dependent on Yalin's defini- 
tions: ''Dunes have a more or less regular two-dimensional shaped! and 
"Ripples have an irregular three-dimensional shape.t' 

As a supplementary means of distinguishing between ripples and 
dunes, velocity profiles fete taken when possible. The intention was 
to compare log-log plots of these profiles (which appear in the appendix) 
with those proposed by Simons and Richardson (ref. 11, p. J26) as being 
typical for particular bed configurations. However, due to the physical 
limitations of the velocity probe, no measurements were taken in the vicin- 
ity of the bed (within approximately 0.03 feet). As such, this supple- 
mentary means of bed form identification was not reliable ale val to be 
abandoned. 

A supplementary means of. identification was also sought to distin- 
guish between dunes and the beginning of the transition to upper regime. 
It was noted during experimentation that as the flow conditions approached 
those at transition, the resistance to flow decreased abruptly, resulting 
in the ponding of water at the foot of the flume, and requiring an adjust- 


ment to the depth to restore uniform flow conditions. Conversely, when 


1 M.S. Yalin, Geometric Properties of Sand Waves, ASCE Journal of the 
Hydraulics Division, Volume 90 HY5, September 1964, p. 115. 


2 ibid. 
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the transition was approached from the upper regime side, an abrupt 


increase in resistance to flow occurred causing a reduced depth at 


the foot of the flume. The occurrence of these abrupt changes in the 


resistance to flow were subsequently employed to aid in identification 


of the bed configuration. 


3-4 Procedure 


A satisfactory laboratory procedure was developed early in the 


testing program and was used throughout most of the experimentation. 


It consisted basically of the following steps: 


2 


Q and S were set at values which would produce a condition 
near transition From one bed configuration to another. 

The flume was filled slowly so as to not disturb the antece- 
dent bed configuration. 

The run was commenced and h was adjusted to obtain uniform 
flow over the bed. 

Depth was readjusted to compensate for any change in the 
resistance to flow due to bed form development. 

Uniform flow conditions were maintained for a sufficient 
length of time to allow complete development of bed forms. 
Photographs were taken of the flume walls indicating the 
shape and dimensions of average bed forms. 

S, Q, C, duration of run (and velocity profiles, if the 
depth was sufficient) were measured. 

The run was terminated, the bed profile plotted and the 
mean depth was determined. 


The bed was slowly drained (so as to avoid disturbing the 
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bed forms) and photographs were taken which indicated 

the shape, size and distribution of bed forms on the 

bed. 
This procedure was followed with various types of antecedent bed configura- 
tions (including plane bed). 

Runs were also made to determine the possibility of plane bed with 
slow (rolling) motion. The procedure employed for these tests was basi- 
cally the same as that previously described. However, the bed was initially 
brought to a plane condition (hydraulically), then screeded to eliminate 
any Pe Uineiinface disturbances, such as ripples. The flume was then 
Operated at a condition approaching initiation of motion. As the bed 
material commenced motion, uniform conditions were maintained and the 
bed configuration observed over as long a period of time as possible. 

Suggestions have been made in the discussion section which require 


slight changes in the apparatus and consequently in the procedure itself. 
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CHAPTER IV 


EXPERIMENTAL RESULTS 


h-] Presentation of Raw and Compiled Data 


Table 4-1 contains the basic observed quantities which were col- 
lected during thirty-one runs. - In addition to these measurements, other 
data were collected: (1) velocity profiles; (2) time variation of sedi- 
ment discharge; (3) time required to reach equilibrium; (4) geometry 
of bed forms. 

The basic data have been reduced to the necessary computed quanti- 


ties and are contained in Table 4-2. . 


h-2 Description of Runs 
h-2-] General 


The experimentation resulted in bed configurations of four basic 
categories: (1) ripples near transition to dunes; (2) dunes near 
transition from ripples; (3) dunes near transition to plane bed; (4) 
plane bed near transition from dunes. No runs were made to determine 


the phase boundaries within upper regime. 


4-2-2 Ripples. 

A total of nine runs were made which resulted in a ripple bed. 
This configuration was developed from both plane bed (screeded by hand) 
and dune bed conditions. The ranges covered by the most important quanti- 
ties were: deena number (0.016 to 0.032), slope (0.00075 to 0.00115), 


h/D (232 to 399) and C (0.2 ppht to 5.3 ppht). Visual identification 
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TABLE 4.2 QUANTITIES COMPUTED FROM EXPERIMENTAL DATA 
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was the primary means of classifying the bed forms as being ripples. 

As such, numerous photographs were taken during and after each run so as 

to obtain a visual record of the bed configuration. The ripple configura- 
tion consisted of forms, usually smaller than dunes, which did not span the 


flume breadth (transversely). This is illustrated by Plates 4-1 to 4-4. 


4-2-3 Dunes. 

Seventeen of the runs made resulted in a duned bed. This configura- 
tion was developed from both ripple and upper regime (plane) Bi ecedent 
conditions. The ranges encountered in the basic quantities were: Froude 
number (0.018 to 0.515), slope (0.00085 to 0.0064), h/D (63 to 373) 
and C (0.5 ppht to 213 ppht). 

The dune configuration was visually unique in that the forms were 
generally larger than ripples and tended to extend across tne entire width 
of the flume. This feature is evident in Plates 4-5 to 4-10. 

Although the longitudinal shape of the dune and ripple forms were 
quite similar, the dunes near transition to upper regime seemed to be 
somewhat more elongated (this becomes evident in a comparison of the pie 


beds of Plates 4-11 to 4-14 with the dune beds of Plates 4-15 to 4-18). 


h-2-h Plane Bed 


Five runs were completed in the plane bed configuration. The 
antecedent conditions included duned and plane beds. The basic quantities 
varied over the following ranges: Froude number (0.128 to 2.212), slope 
(0.0066 to 0.0075), h/D (52 to 155) and C (226 ppht to 1063 ppht). 

The plane bed was easily distinguished, visually, as is evident 


in Plates 4-19 to 4-21. In some instances, however, it was discovered 
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PLATE 4-1. RIPPLE BED OF RUN 
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PLATE 4-3. RIPPLE; BED; OF. RUN 


PLATE 4-4. RIPPLE BED OF RUN 
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PLATE 4-7. DUNE CONFIGURATION OF RUN 


PLATE 4-8. DUNE BED OF RUN 22 
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PLATE 4-9. DUNE CONFIGURATION OF RUN 22 


PLATE 4-10. DUNE BED OF RUN 24 
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PLATE 4-13. RIPPLE FORMS OF RUN 


PLATE 4-14. RIPPLE FORMS OF RUN 
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PLATE 4-15. FORM OF DUNES NEAR TRANSITION FROM RIPPLES (RUN I!) 


PLATE 4-I6. FORM OF DUNES NEAR TRANSITION FROM RIPPLES (RUN |6) 
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PLATE 4-17. FORM OF DUNES NEAR TRANSITION TO PLANE BED (RUN 22) 


PLATE 4-18. FORM OF DUNES NEAR TRANSITION TO PLANE BED (RUN 24) 
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PLATE 4-19. PLANE BED OF RUN 20 


PLATE 4-20. PLANE CONFIGURATION OF RUN 19 
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PLATE 4-21. PLANE BED OF RUN 19 
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h3 
upon completion of the run that some patches of very low, flat dunes existed 
near the flume walls (see Plates 4-22 to 4-24). This was generally encoun- 


tered when any tendency towards meandering within the flume walls occurred. 


4-2-5 Plane Bed at Low Charge 


Although several runs were made in an attempt to obtain plane bed 
with rolling motion (at low charge) this configuration could not be main- 
tained. In all the runs the bed, which had initially been screeded to a 
plane condition, developed a ripple formation. An example of this phenomenon 


is shown in Plate 4-25. 


4-3 Supplementary Data 


During the ete of experimentation certain data were collected 
which were of interest but not of direct importance to this study. A plot 
of the time variation of sediment discharge (as determined by periodic 
sampling) is included in the appendix as Figure A-]. Profiles of the 
bed at various elapsed times after commencement of a run (to determine 
the time required to reach equilibrium) are illustrated by Figures A-2 to 
A-4 respectively, in the appendix. Finally, Figures A-5 to A-16 represent 
the previously discussed velocity profiles, and Figures A-17 to A-27 are 


the bed profiles which were plotted at the termination of each run. 


h-h Presentation of Results 


The basic quantities, as derived in Chapter V, are graphically 
illustrated by the two-dimensional plots of Figures 4-1 to 4-7. In 
addition to the experimental points, the transitional boundaries as pro- 


posed by Cooper (3) are also plotted. 
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PLATE 4-23. SHALLOW DUNE FORMS NEAR WALL (RUN 27) 


PLATE 4 -24. SHALLOW DUNE FORMS NEAR WALL (RUN 28) 
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PLATE 4-25. DEVELOPMENT OF 
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CHAPTER V 


ANALYSIS AND DISCUSSION 


5-1 Prediction of Bed Forms and Their Behavior 
5-1-1 General 


As stated in the introduction, it is desirable to establish a 
mathematical model which, when provided with a minimum of independent 
variables, would predict bed configuration. However, before such a 
model can be constructed, all the variables which are related to the 
phenomenon must be specified. Thus, the following physical properties 
must be present in the statement of the problem (for a flume study): 

(1) flume properities (breadth, b). 


(2) fluid properties (density Pes and viscosity, v). 


(3) coefficient of gravitational acceleration, g. 


(4) bed material properties (median diameter, density, Ph» 


D 
50% 
and a factor to account for particle shape, gradation, etc., 


X). 


(5) the flow properties (discharge, Q, mean velocity, V, slope, S, 


depth, h, and rate of sediment transport, C). 
(6) bed form properties (type, wavelength, A, wave height, A). 
Therefore, the general problem statement is: 


VesceG. As 0 -= Fn (Q, h, Pee v, b, Pb» D X, g) (5-1) 


50’ 


Thus, a dimensionally complete statement must contain ten variables. 


According to dimensional analysis theory (ref. 2), the presence of three 
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basics (mass, length, and time) fixes the number of dimensionless parameters 
at seven. If the starting variables chosen are: 


fn (0) » Pr» X, b, NyaveeD 0? Gy V, C) = 0 (5-2) 


>) 


then one possible set of parameters is: 


fn, WCE! Toth 9, oC Alvient Romrex fab) of= 5<0 53) 
50 Pr p 
y2 
Where FR' = PF and V. = 3 /vg D9 


a(,-pe)h sei at 


The last four parameters in the list are assumed to be either constant or 
ineffectual (although V, would vary with any change in bed material, it 
remained constant during this study). Thus, a simple, practical statement 


which describes the mathematical model is: 


treaties. -C),,-s40 7 (5-4) 


Similarly, if the starting variables are: 
fn (po, Peo X, b, h, v, D, g, S, C) = 0 (5-8) 


then a second statement which defines the system is: 
Frese cy. h/D). =O (5-6) 


Thus, the mathematical model may be represented graphically with a three- 
dimensional co-ordinate system upon which the bed configuration boundaries 


may be indicated. 


5-1-2 Comparison of Prediction Schemes 


In order to verify the mathematical model thus obtained, it would 


be useful to be able to compare it with other models derived from previous 
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analyses. Unfortunately, the most widely accepted prediction scheme (that 
of Simons and Richardson, as shown in Chapter I1) employs a unique measured 
quantity, fall diameter, and is not dimensionless. Similarly, Liu's model 
(see Chapter I!), although dimensionless, employs the variable w, fall 
velocity, which inhibits direct comparison of solution curves. As such, 
only the Simons and Richardson criterion is useful as a basis for compari- 
son with the experimental results. 

Figure 5-1] illustrates the experimental data, as well] as curves 
which best fit the Gilbert data, in relation to the Simons and Richardson 
curves. By inspection, it becomes obvious that the data and curves are 
in agreement only with respect to the sequence of bed form development. 
Admittedly, the Gilbert data contain no points which are classified as 
ripples; however, it appears as though the transitions in bed configuration 
occur at significantly lower values of stream power (t VY) than those pre- 
dicted by Simons and Richardson. 

In order to gain a better understanding of the Simons and Richardson 
criterion for prediction of bed configuration, it becomes necessary to 
examine its derivation from a dimensional analysis point of view. The 
original problem statement in the derivation (1964) is: 

PU Tn 195) D5 ly) Wy. Py Oe (5-7) 
The last two variables may be eliminated from the expression by specifying 
a range of experimental conditions. Thus, the expression becomes: 

Formay=o arn. (S,1-D;- hye} (5-8) 
Simons and Richardson state that ''The variables were not grouped into 


dimensionless parameters for discussion purposes to avoid masking the 
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essential role of any one of eHemetes and derived the final solution: 


Forme = eine (t:V, D 
Oo 


fall’ Ck 


The Simons and Richardson solution would contain the proper number 
of parameters, if the quantities involved were dimensionless. However, by 
lumping several quantities together with V to form a physically meaningful 
parameter, TOY the authors appear to have eliminated these quantities 
from further consideration in the problem. This is analagous to attempt- 
ing to reduce a statement of the form: 

X = fn (Y, Z, W) 7 (5-10) 
(which results in a solution surface in three dimensions, to the form: 

XYZ = fn (W) (Sait) 


(which may be represented by a two-dimensional co-ordinate system). 


5-2 Bed Form Dimensions 


lf the list of starting variables from Section 5-1-1] were as 


follows: 
Deno Ga vs = 0 meats 
Fn (o,, Peo Neat saitle sy 50 g, V r) (5 ] ) 
another set of parameters may be derived: 
ame AE VAN VAS) eat ae (5-13) 


Although the analysis of wave length and wave height is beyond 
the scope of this study, a special effort was made to collect data which 
might be useful for such an investigation at a later date. As such, 


the average height and length of the waves for each run were plotted 


1. Simons, D.B., Richardson, E.V., and Nordin, C.F.JR., "Sedimentary 
Structures Generated by Flow in Alluvial Channels'', The Society of 
Economic Paleontologists and Mineralogists, The American Association 
of Petroleum Geologists, Special Publication No. 12, p. 47. 
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against FR' as shown in Figures 4-6 and 4-7. 

As the bed form geometry data were copettered to be unreliable, no 
direct comparisons could be made with the work of other authors. It is 
significant, however, that Yalin (13) considered only the parameters h/D 
and A/D in his plot which predicts bed form configuration (see Figure 


5-2). 


5-3 Discussion 


Due to the nature of this experimental study, a number of problems 
were encountered. Basically, the difficulties arose from the large number 
of variables involved in the phenomenon of sediment transport. It was 
one thing to assume that a transition in bed configuration might occur 
near a point with certain flow conditions, and yet another thing to attempt 
to stabilize the experimental system at that point. Although only two 
quantities could be imposed on the system (eg. Q and h), the time required 
for the remaining variables to reach their equilibrium values changed 
from run to run. 

This fact led to some problems related to the actual operation of. 
the flume. The system was allowed more than adequate time (up to forty- 
eight hours) in which to reach equilibrium. This resulted in a substantial 
amount of evaporation which affected the depth of flow (by reducing the 
volume of water in the closed system). Thus, one of the supposedly 
"imposed'' quantities was actually slowly changing throughout the run. 

The amount of time required to complete a run was also affected 
by rapid adjustments in Q, S or h. An increase in discharge or a decrease 
in ieee or depth caused the downstream end of the bed to momentarily 


assume control over the flow. This invariably caused scour for some 
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distance upstream and eventually resulted in deposition of the removed 
material at the upstream end of the flume. Thus, more time was required 
to reach the desired flow conditions, and to restore the bed to its 
original condition (approximately parallel to the flume bed). 

This characteristic instability of the sand bed in the vicinity 
of the inlet to the pump (at the foot of the flume) also affected the 
total load samples. The time variation of total sediment discharge, as 
shown in Figure A-] in the appendix, resulted in large deviations from 
the mean. As these surges could not be controlled, it was decided to 
accept them as being inevitable and take as large a number of samples 
as possible. 

As mentioned in Section 5-2, the bed form geometry data obtained 
were unsatisfactory. This resulted from the bed plotter's lack of sensiti- 
vity which often led it to smooth out the extremities of bed forms and 
even caused it to skip over some of the smaller forms. Furthermore, 
it was impossible to use the plotter during the run (which would have 
greatly simplified the determination of depth) as the relatively large 
head on the probe caused scouring of the bed and resulted in a false 
profile. This inability to obtain reliable information about bed form 
geometry was especially critical. Due to the size of bed material 
employed and to the range of h/D encountered, the dimensions of ripples 
and dunes were of the same order of magnitude (which hampered bed form 


identification). 
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CHAPTER VI 


CONCLUSIONS AND RECOMMENDATIONS 


6-1] Summary and Conclusions 


The object of this study has been to define, as precisely as 
possible, the conditions at which transition in bed configuration occur. 
A review of existing knowledge indicated that many theories are limited 
by their analytical approach and that others have not considered all of 
the variables relevant’ to the sediment transport phenomenon. As such, 
an attempt has been made to develop a mathematical model and establish 
its validity by.comparison with experimentally derived points. 

The mathematical model was developed by means of dimensional 
analysis so as to ensure the completeness of the final result and to 
eliminate as many simplifying assumptions as possible. Its final form 
may be expressed as: 

Th AG re Ae (6-1) 
and cree 950C, 0/0) a= 7 0 (6-2) - 
The model, upon which the phase boundaries proposed by Cooper (3) were 
plotted, was evaluated by placing the experimental points on two-dimensional 
plots of the parameters stated in Equation 6-1. The boundaries, which were 
derived from a voluminous amount of flume data, were, for the most part, 
substantiated by the experimental data from this study. 

However, it is felt that the Cooper curves at the transition from 
dunes to upper regime should, in fact, be a zone. One significant result 
of this change would be that two distinct families of C contours would 


exist (on the Fr' versus h/D plot, say) in the transition zone. Figure 
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6-1 shows some hypothetical contours of C superimposed on the Fr' versus 
h/D plane. In the transition region, C is expected to be greater if the 
antecedent condition in the bed form development were plane bed rather 
than dunes. 

Furthermore, it would seem from the experimental points that the 
transition from ripples to dunes occurs at lower values of FR! than those 
proposed by Cooper. This boundary could perhaps be adjusted as indicated 
by Figure 6-1. 

The lack of a plane bed configuration at low charge for sand size 
material indicates that the bed material properties which were eliminated 
by Cooper may play an important role in determining some phase boundaries. 

The difficulty encountered in distinguishing between ripples and 
dunes may Poaicere that, in fact, there is no distinction between the two 
configurations in the experimental range of h/D. However, as there is also 
an imperceptible difference in resistance to flow and sediment transport, 
perhaps there is no need to distinguish between these bed forms in the 
experimental range of h/D. If such a distinction were necessary, an adjust- 
ment of the Cooper ripple-dune boundary (as indicated in Figure 6-2) might 


be useful. 


6-2 Recommendations 


As this has been an experimental study, the recommendations deal 
mostly with proposed modifications to the apparatus. These changes are 
suggested in light of the experience gained in operating the existing 
system: 

(1) Cooling coils could be added to the return line to reduce 

evaporation or an automatic system could be developed 


(perhaps by means of a float arrangement) to control the 
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quantity of water in the flume. 

(2) The flume bed and walls could be extended beyond the inlet 
to the sand pump in order to provide a damping effect to 
any sudden adjustments in Q, S or h. 

(3) A better means should be found for obtaining bed form dimen- 
sions (an optical bed plotter is currently under development 
at the University of Alberta). 

(4) A pitot tube could be used to obtain more precise velocity 
profiles over the various bed configurations. This system 
would necessitate slow backflow to prevent its being plugged 
with sediment particles. 

With regard to further experimental work, it is becoimendet that 
more efforts, Peecit ion to the improved apparatus mentioned above, be 
applied to the determination of the relationship between bed form geometry 
and flow conditions. It is also recommended that the significance of 
bed material properties in determining the location of the plane bed (low 
charge) to ripple or dune phase boundaries be investigated. 

Furthermore, in order to substantiate the existence of a zone of 
transition between dunes and upper regime, it is suggested that further 


testing be done with C constant and maintaining a uniform bed configuration. 
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